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ABSTRACT 



A fast, accurate and portable data reduction system 
was developed for the XR-3 Captured Air Bubble testcraft 
being evaluated at the Naval Postgraduate School, Monterey, 
California. The system consists of four units. A magnetic 
tape recorder is used for data play back. A signal condi- 
tioner unit with a built-in analog- to-digital converter was 
developed and is used to filter, amplify, sum and further 
prepare the data for transmission to either a strip chart 
recorder or a digital X-Y plotter through a Monroe 1880 
calcu la tor . 

Preliminary use of curve fitting techniques are dis- 
cussed; ^calculator programming and the various problems 
and solutions encountered in the development of the system 
are described. 
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I 



INTRODUCTION 



Associated with most experimental research projects 
is the requirement to obtain, evaluate, Interpret, and 
analyze large quantities of data in order to draw meaning 
ful conclusions. A manual data reduction method is both 

I 

time consuming and prone to error. On the other hand, an 
automatic data reduction system is more accurate and can 
significantly shorten the man-hours required to reduce 
raw experimental data into a form which can be quickly 
utilized by the project engineers. 

An automatic data reduction system was needed for the 
XR-3 Captured Air Bubble (CAB) testcraft being evaluated 
by the Naval Postgraduate School (NPS), Monterey, Califor 
nia. This craft is being operated at Lake San Antonio, 
approximately 100 miles from the school. Test data are 
recorded on a 14-track Pemco, Model 120B,. magnetic tape 
recorder. The method used to reduce the data has con- 
sisted of transporting the tape unit back to NPS, using 
a digital voltmeter and/or a strip chart recorder to ob- 



tain a visual data 


presentation , 


and then 


making 


hand 


performance plots . 


The ratio of 


data acqu 


isitio 


n time to 


reduction time was 


about 1 to 5. 


The need 


for a 


faster 



system arose from the decision to acquire a facility that 
would permit the test operators to remain overnight at th 
test site. In order to maximize the efficiency of 
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operating on-site over a several-day period, the results 
of each day*s testing must be reduced and analyzed before 
a schedule of tests can be prepared for subsequent test 
runs* Therefore, the time to reduce the data must be 
shortened significantly. 

The automatic data reduction system described in this 
paper fulfills these objectives. It provides for a better 
utilization of operating time at the test site by provid- 
ing a quick, accurate and portable means of evaluating the 
performance of the testcraft. 
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BACKGROUND 



A. NATURE OF THE PROBLEM 

Both the data acquisition system and the data reduc- 
tion methods for the XR-3 have been constantly improved 
since the testcraft arrived at the Postgraduate School in 
March, 1970. Initially, performance data were limited to 
thrust, velocity and plenum cavity pressure. These para- 
meters were hand recorded and hand reduced to produce 
rough studies of thrust versus velocity and optimum plenum 
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engineering units, such as knots for velocity vice the 
millivolts (mV) output from the tape recorder. The system 
had to be portable, in order to be moved in and out of the 
data reduction facility, and also small and capable of be- 
ing operated in a limited space. The equipment had to be 
simple and easy to use, requiring few inputs from the 
operator. One of the major considerations was that the 
system had to be expandable, which would allow for the 
implementation of future requirements. 

The presentation of the results was considered an 
important factor. Three-digit data accuracy was sufficient 
both for plotting and for the printed output since the 
raw data are accurate only to this order of magnitude. The 
data would require being scaled twice, once for the plot, 
changing millivolts to inches on the graph, and once for 
printing, if desired, changing millivolts to physical units. 
The system also had to be capable of performing operations 
such as summing inputs from two or more channels of the 
tape recorder. One such application is in determing the 
total drag on the testcraft at a constant speed. In this 
condition the drag equals the thrust, which in turn equals 
the sum of the thrust from both the port and the starboard 
engines . 

All of the phenomena being tested occur at a fairly 
low frequency, no greater than about 0.5 Hz. For this 
reason the sampling rate need not be much greater than 
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one Hz according to Shannon's sampling theorem (also known 
as Kirchhoff's sampling theorem). However^ to allow for 
future expansion a faster sampling rate would be needed. 

B. MEETING THE OBJECTIVES 

A Monroe model 1880 scientific programmable printing 
calculator and a Monroe model PL4 digital plotter were 
chosen to be the basic components of the new data reduction 
system (see Figure 1). Because of its small size and rela- 
tively large 'instruction set^ the calculator seemed best 
suited to perform any manipulations necessary to prepare 
the data for the plotter and also to provide a paper tape 
printout of the data converted to the proper physical units. 

An analog- to-digital (A/D) converter was necessary to 
transform the analog tape recorder signal to a digital form 
compatible to calculator input. A contract was let to 
Santa Cruz Engineering of Santa Cruz^ California in April, 
1973 to design and build an interface which would make it 
possible for the calculator to select a tape recorder 
channel, then convert the signal from that channel to a 
digital input to the calculator. Calculator software (pro- 
gram) would convert this input into the various forms 
necessary to accomplish the remaining tasks. A delivery 
date of June 1973 was set; however, an extension was re- 
quired until late July 1973. 

As an intermediate method of data reduction while 
waiting for the interface, various methods of curve fitting 
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Figure 1. Monroe 1880 Calculator and PL4 Plotter 
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III. CURVE FITTING 



A. GENERAL METHODS 

There are various methods available to engineers for 
fitting data to an algebraic relation that lends itself 
to easy manipulation and plotting and that can have a 
physical significance. The methods explored for applica- 
tion to the data reduction system were the following: 

1. Piecewise - polynomial interpolation method 

2. Difference method 

3. Least-squares approximation by orthogonal poly- 
nomials method 

4. Least-squares regression method 

1 . Piecewise -Polynomial Interpolation Method 

The cubic spline method of piecewise-polynomial 
interpolation was not satisfactory for use by either the 
calculator or the IBM 360 computer for two reasons. 

First, the method required knowledge of the slope of the 
curve at the initial and final points. This information 
is not known and only a gross estimate could be achieved. 
Secondly, a spline fit requires the curve to pass through 
each of the data points. Recording error and signal 
noise are not smoothed by this method and as a result an 
inaccurate model of the physical event results. A spline 
fit is best suited to determine polynomial approximations 
for exact data whose slope has discontinuities [Ref. 2]. 
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Difference Method 



The difference method likewise does not lend it- 
self to computer application; however, it is an excellent 
tool for determining the degree of fit required for a 
least-squares type method. When the difference method is 
employed, the data must be gathered such that the depen- 
dent variable is obtained at equal increments of the in- 
dependent variable (X) . A difference table is formed 
taking iterations until an essentially constant difference 
is reached throughout a given iteration. The order of the 
polynomial to give the best fit is the same as the number 
of iterations to reach this constant difference [Ref. 3]. 

That is, if the second iterations produce a constant value, 

2 

then the model will be of the form y = a^ + ^ - ^ 2 ^ 

The coefficients of the polynomials are solved for by 
using the Gr egory-Newton interpolation formula: 

r(r-l)(A^y ) r (r-1) (r-2 ) ( A^y ) 
y = + r(AyJ + ^ 

r (r-1 ) , , (r-n) ( A^y ) 

+ + 2 _ 

where 

X - X 

o 

r = 

w 

w is the increment between successive values of the 
independent variable x, and 

y , x are the initial values of the variables 
o o 
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EXAMPLE 1 DIFFERENCE METHOD 



As an example of this method, the data for a third 
order equation y = 5 + 2x^ are presented in Table I. 
Since the constant difference occurs in the third itera- 
tion, the desired fit will be given by a third order 
polynomial. In addition, since the third iteration re- 
sulted in a difference value with zero variance, the fit 
will be exact. 



TABLE I 

Successive Differences for Difference Method of Curve 
Fitting . 



Given Data 
X y 


Successive Differences 


1st 

Iteration 

Ay 


2nd 

Iteration 


3rd 

Iteration 


4 th 

Iteration 




0 5 

2 21 

4 133 

6 437 

8 1029 

10 2005 


Ay^= 16 
Ay^= 112 
Ay^= 304 
Ay^= 592 
Ay^= 976 


A^y = 96 

2 

A y^= 192 
288 

2 ^ 

A y^= 384 


A ^ y = 9 6 

3 ^ 

A y = 96 
A^y^= 96 


A y^= 0 



^0= ^0= 5, W = 2 



r = 



x-0 



X 

2 



and 



y = 5 + f(16)+ |(| -1) (|^)+ |(| -1) (| -2) (|i) + 0 



y = 5 + 2x' 
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which is the exact equation from which the data were de- 
rived. To determine the order of a model from data which 
is not exact, each column of the difference table must be 
inspected for an essentially equal distribution of plus 
and minus signs. The number of iterations to that point 
determines the order of the resulting polynomial fit. 

The remaining procedure is as discussed in the preceding 
example . 

3 . Least-Squares Approximation by Orthogonal 
Polynomials 

The. method of least-squares using orthogonal poly 
nomials could not be used as a calculator solution to the 
problem since more storage is required by the method than 
is available in the calculator. However, it does provide 
a very efficient and accurate solution when used on the 
IBM 360/67 computer. 

The method consists of solving for successively 
higher order orthogonal polynomials, such that the sum of 
which satisfies the requirement that 
N 

[f(x ) - p(x )]^ where 

, n n 



f(x) = 


exact 


solution 


p (x) = 


polyn 


omi al approximation 


N = 


order 


of the polynomial 



is as small as possible [Ref. 2]. A computer program was 
developed for use on the IBM 360/67 in which the X and y 
data points, the number of points and the degree were 



\ 
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entered. The output consisted of two parts: The first 

contained : 

1) The coefficients of the fitting polynomial 

2) An estimate of the error of each coefficient 

3) Sum of the square of the deviations 

4) The F-ratio 

5) The Chi Square goodness of fit statistic 

6) Degrees of Freedom of the data, and 

7) A list of the data points, the evaluated poly- 
nomial at each value of X, and the differences 
between the two y values. 

Part two is a plot of: 

1) The data points, 

2) the fitted curve, and 

3) the normalized error. 

The computer program and sample output can be found in 
Appendix B. 

4 . Least-Squares Regression Method 

The least-squares regression method of data fitting 
is the most familiar and the most widely used by engineers. 
It is used by hypothesizing the existence of a polynomial 
that will best fit the data as in the case of orthogonal 
polynomials. The degree of the fitting polynomial must be 
determined by another method or by analysis of the results 
in order to determine the "best” fit. In the case of an 
Nth order fit, the fitting polynomial has the form: 

2 3 N 

y = a + bx + cx + dx + + nx 
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Additional equations are needed to solve for the coeffi- 



cient a. 


b/ c^ • t 


n. These 


are ge: 


plying the equation by 


2 3 

X, X , X , 


etc. , 


variable 


on both sides 


to obtain: 




Ey 


= Ma +b X 


+ C 


+ d E : 


E*y 


= a ^ X +by^ x^ 


y* 3 

+ C JL^ X 






= a 


V 4 
+ c X 


-0 E^ 



n+1 



n+2 



( 1 ) 



/ y = ^ X +b Z-r X +c Z-^ X +d Z-r x +...+n Z-^; 



2n 



where 



M is the number of points. 

In matrix notation Equation (1) can be written: 
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EXAMPLE 2 Leas t“squares regression method 

As an example, the data of example one will be fitted 
to a third order polynomial. The data in Table I can be 
used to generate the data in Table II. 



TABLE II 

Least-Squares Data 







2 




3 


2 


4 


3 


5 


6 


X 


y 


X 


xy 


X 


X y 


X 


X y 


X 


X 


0 


5 


0 


0 


0 


0 


0 


0 


0 


0 


2 


21 


4 


42 


8 


84 ■ 


16 


168 


32 


64* 


4 


133 


16 


532 


64 


2128 


256 


8512 


1024 


4096 


6 


437 


36 


2622 


216 


15732 


1296 


94392 


7776 


46656 


8 


1029 


64 


8232 


512 


65856 


4096 


526848 


32768 


262144 


10 


2005 


100 


20050 


1000 


200500 


10000 


2005000 


10000 


100000 


Z 1 30 


3630 


220 


31478 


1800 


284300 


15664 


2634920 


51600 


412960 



Substituting the data from Table II into the Equation (1) 



3630 = 6a + 30b + 220c + ISOOd 

31478 = 30a + 220b + 1800c + 15664d 

284300 = 220a + 1800b + 15664c + 51600d 

2634920 = 1800a + 15664b + 51600c + 412960d 



(4) 



Solving (4) for the coefficients gives 
a = 5, b = 0, c = 0, d = 2; 
therefore the resulting equation is 

. ^ 3 

y = 5 + 2x 

which is the exact solution to the data being fitted. 
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B. CALCULATOR APPLICATIONS 

Experimentation with representative data on the IBM 
360/67 computer showed that a polynomial of seventh order 
was required to best-fit a curve smoothly through the 
data points. This was also verified using the difference 
method. The least-squares regression method of curve fit- 
ting was implemented on the calculator solving for a 
seventh order polynomial. Various problems were en- 
countered while trying to solve Equation (1) for the un- 
known coefficients (a,b,c, •••/ n) because the A matrix 

was highly ill-conditioned. Two methods were used to try 
to solve the equation. The first was a direct approach 
using triangular factorization with forward and back sub- 
stitution. The second was an iterative method using a 
Gauss-Seidel iterative scheme with and without successive 
over-relaxation (SOR) . Examples of these methods can be 
found in reference 4. 

The triangular factorization did not produce accept- 
able results because of the loss of significant digits. 

The calculator solution differed so extremely from the 
computer solution and the plot of the resulting polynomial 
equation using the calculator solution bore so little 
resemblance to actual performance plot that other methods 
had to be tried. Scaling of the input data had no correct 
ing effect on the results using this procedure. 

The Gauss-Seidel iterative method occasionally would 
provide the proper coefficients for the fitting polynomial 
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however, this was also an unsatisfactory method as a 
general solution to the problem. See Figure 2 for the 
flow chart of the calculator program implementing the 
Gauss-Seidel iterative method. The initial guess that 
was used to arrive at the proper solution was extremely 
critical. An improper guess would lead to instability 
and the coefficients wo'uld diverge from the proper solu- 
tion. Even when stable and converging, the time required 
for the calculator to reach a solution was excessive. 

An iteration required about 20 seconds and, depending on 
the initial guess, a problem could require over 200 itera- 
tions to arrive at a solution which would closely approxi- 
mate the physical phenomenon. The introduction of a 
relaxation factor helped; however, like the initial guess, 
the choice was critical and no way was developed to aid 
in the choice of either an initial guess or a relaxation 
factor . 

The time problem and the initial guess and relaxation 
factor problems caused this method to be unacceptable as 
a general solution to the problem of curve fitting on the 
Monroe 1880 calculator. 
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Figure 2. 



Flow Chart of Calculator Curve Fitting Routine. 
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Figure 2. (continued) 
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IV. DATA REDUCTION SYSTEM 
A. SYSTEM COMPONENTS 

The data reduction system is composed of the fol 
five major components: 

1. Signal selector and conditioner unit 

2. A/D converter and calculator interface modul 

3. Monroe 1880 calculator 

4. Monroe PL4 digital plotter 

5. Hewlett Packard model 7100B strip chart reco 

The data in the form of an analog signal are pla 

back from all 14 channels of the Pemco magnetic tape 
corder. These signals go into the signal selector a 
conditioner unit from which the operator chooses the 
that he wants each of the signals to follow. He can 
up to eight signals to the A/D converter module also 
located in the signal selector and conditioner unit, 
in turn feeds digital raw data to the Monroe 1880 ca 
tor for manipulation and plotting, or he can send th 
analog data to the multi-channel strip chart recorde 
This flexibility is essential to the operation of an 
efficient data reduction system. 

1 . Signal Selector and Conditioner Unit 

The signal selector and conditioner module w 
signed to allow for operator ease in selecting, moni 
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and routing signals from the tape recorder to the other 
components of the data reduction system. All channels of 
the recorder feed into the module through the first 14 
terminals of a CO-ORD Switch 10 x 15 pin, 2-deck Matrix 
Program Board. By selecting the desired output channel, 
the operator can filter, amplify and/or sum signals using 
any of nine filter/amplifier circuits and a summing cir- 
cuit on the 10 outputs of the matrix (see Figure 3). The 
signal voltage can be monitored at any of 24 points using 
a built-in Datel model DM 2000AR Digital Panel Meter and 
a 24 position rotary selector switch. The monitoring 
points include each of the inputs from the 14 tape re- 
corder channels, each of the outputs from the nine filter/ 
amplifier circuits and the output from the summing circuit. 

Signal conditioning was needed to filter out the 
extraneous high frequency noise from the tape recorder 
signals which originated from the XR-3 onboard sensors. 

The noise results principally from vibrations of the test- 
craft during normal operations; consequently each of the 
signals has about the same characteristic frequency noise 
as shown in Figure 4. The filters used are a Butterworth 
low pass type [Ref. 5] with a cutoff frequency of 0.5 Hz. 
The electronic circuit diagram and the signal to frequency 
response curve are shown in Figure 5. The circuit was 
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Figure 3. Signal Selector and Conditioner Unit. 



29 




Figure 4 • 



Sensor Output Signals. 
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Butterworth Filter Circuit 
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maximum signal from the tape recorder must be amplified 
to a five-volt maximum signal for use with the A/D con- 
verter. This is to ensure the best possible data 
accuracy when using the A/D converter and calculator in- 
terface module. Since each of the sensor/amplifier/ 
signal conditioner circuits on the testcraft does not 
maintain a null output for a zero signal input, a zero 
offset voltage circuit was included. Inasmuch as a full 
scale input does not always have an exactly 1.0 volt out- 
put, a variable gain resistor was placed in the feedback 
circuit of the amplifier to allow the operator to adjust 
the gain to have a 5.0 volt output for a full scale input. 
See Figure 6 for the amplifier circuit. 

Two of the filter/amplifier circuits are added 
then halved through a summing circuit using another Burr 
Brown 3440J amplifier (see Figure 7). This circuit is 
used in obtaining total thrust from the individual port 
and starboard thrust components. 

Electrical power is supplied to the Burr Brown and 
Philbrick amplifiers from a, Philbrick/Nexus NPS-300 +15v 
-15v power supply. Zero offset voltage is supplied by a 



SRC MOD 


3564 powe 


r 


supply and the 


Datel 


Digi 


tal Panel 


Meter is 


powered 


by 


a Acopian MOD 


10J75 


lOv 


dc power supply 


reduced 


to 5v dc 


by 


a National Sem 


icond 


uc tor 


LM340K-7805 



5v dc regulator. 
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Figure 6. Amplifier Circuit. 
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2 • Analog-to Digital Converter and Calculator 
Interface Module 

The analog- to-digital converter and calculator in- 
terface module, also called the A/D interface, was pur- 
chased from Santa Cruz Engineering in late July 1973o Many 
modifications were required to permit it to be used with 
the Monroe calculator and the other components of the data 
reduction systemo These modifications are discussed in 
Section V, i 

The A/D interface is selected through connections 
on the front of the signal selector unit. As many as eight 
signals can be connected and generally the outputs from the 
filter/amplifiers are used as the inputs. The calculator 
program selects which input signal is to be converted by 
sending a three-digit octal code from the calculator IX 
register to a TTL/MST 74164, S-bit, serial to parallel con- 
verter in the interface. The first three bits are then 
sent to a RCA COS/MOS analog multiplexer CD4051A for the 
channel selection. The multiplexer receives the input 
signals from all eight input channelso It acts as a binary 
controlled analog switch by selecting which channel is to 
be outputted to the A/D converter. The A/D converter is 
then triggered by a TTL/MONOSTABLE 9602 multivibrator. 

The A/D converter is a high speed, highly accurate 
and stable Analogic MP 2410 unit housed in a small 4 inch 
X 2 inch x 0.4 inch packageo It is capable of ten-bit 
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resolution; however, only eight bits are used for com- 
patibility with the Monroe 1880. It is capable of ranges 
of -lOV to +10V, -5V to +5V and OV to +5V with output in 
one’s complement, two's complement, offset binary or 



unipolar binary. The OV 


to 


+5V range using 


unipol 


ar 


binary was selected 


for 


the 


data reduction 


system 


appli- 


cation. The output 


from 


the 


A/D converter 


then go 


es into 



a TTL/MSI 7491, eight-bit, serial-in serial-out, shift 
register. Then upon another command from the calculator, 
the binary signal is transferred from the shift register 
in the interface to the IX register of the calculator for 
further processing. 

Timing for the interface is provided by the A/D 
converter for the internal function and by the calculator 
for data transfer. Power is provided by an Analogic 
AN3001 modular power supply outputting +15V, -15V and 

+5V dc. Timing and transfer logic is provided by FAIR- 
CHILD TTL/LSI 7400 Quad two-input NAND GATE and 7410 
triple three-input NAND GATE integrated circuits. The 
logic diagram for the A/D interface can be found in 
Figure 8. 

3 . Monroe 1880 Calculator 

The Monroe Model 1880 Scientific Programmable 
Printing Calculator was chosen as the heart of the data 
reduction system because it combines the features of a 
simple-to-oper a te but powerful scientific calculator with 
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Figure 8. A/D Converter Logic Diagram 



an extensive problem-solving capability. The model being 
used provides 1024 memory locations for program steps (in- 
structions), 128 main data registers and 10 temporary 
(scratch pad) registers (both for data storage) . The main 
data registers can be programmed to hold an additional 
1024 instruction steps instead of data if a larger program 
is necessary. Data is stored to an accuracy of 13 signi- 
ficant digits. 

The operations of the calculator fall under six 
general categories: input, output, storage, control, 

execution and arithmetic computation. All of these opera- 
tions can be programmed from the keyboard either through 
individual function keys or through an ENTER CODE key 
which allows the three-digit octal code of the instruction 
to be entered into memory (see Figure 9). 

A program can also be stored using a magnetic card 
device which is built into the calculator and enables pro- 
grams to be both read into the calculator and written on 
to small magnetic cards each capable of holding 512 program 
steps. Data from the main data registers can also be 
written on the cards and conversely data from the cards 
can be read into the main data registers. Data may be in- 
putted to and outputted from peripheral devices through 
the IX register and an external 36-pin connector to the 
peripheral. Programs can be printed on paper tape with 
the listing containing the instruction number, the three- 
digit octal code and a symbol to identify the instruction 
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Figure 9. Monroe 1880 Keyboard. 
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if one is applicable. Program steps can also be inserted, 
deleted or changed without having to re-enter the entire 
code. This feature is a tremendous aid in correcting 
errors in a program. 

Control functions include such items as sitting 
flags, printing data, and setting the position of the 
decimal point. Program execution is extremely easy with 
the Monroe 1880. Symbolic program addressing makes it 
possible to locate a program anywhere in program memory. 

It permits branching or jumping to any step in a program 



or to the 


beginning of a 


subroutine without 


h 


aving to 


know the 


absolute storage 


location address. 




Branching re 


turns to 


the point of the 


branch after the 


ex 


ecution of a 


reserve ; 


a jump does not. 


Decision making 


is 


made simple 



with the use of a sense switch, a programmable and a key- 
board flag and through the test of a positive, negative 
or zero value of the contents of the main entry register. 
Tests of these devices are used to control branches or 
jumps to other parts of the loaded program. 

The arithmetic functions performed by the calcula- 
tor include such standard keyboard items as add, subtract, 
multiply, divide, square root and reciprocal. In addition, 
other keyboard items such as sine/cosine, arc sine/arc 
cosine, logarithm, statistical summation and factorial are 
included. Many more functions and controls are available 
through enterable instruction codes for which keys are not 
available. These include most of the codes for transferring 
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information to peripheral devices, program control, regis- 
ter manipulation and other arithmetic functions such as 
absolute value and arc tangent. A complete set of instruc- 
tion codes and operating instructions can be found in 
Refs. 6, 7, and 8. 

4 . Monroe PL4 Digital Plotter 

The Monroe PL4 Digital Plotter is a peripheral de- 
vice used in association with the Monroe 1800 series calcu- 
lators. Its purpose is to display graphically data which 
is passed to it as a result of program execution. The graph 
origin may be set anywhere on the plotting surface through 
the use of a calculator scaling routine which also formats 
the data for transfer to the plotter. The plot subroutine 
discussed in Section IVB and shown in Appendix D, determines 
whether the plotter is "busy" (actually plotting) or is 
ready to receive data. If busy, the program loops until 
the plotter is ready at which time it passes the pen command 
and the data for plotting. 

Instructions and data are passed to the plotter 
through the calculator IX register. The four commands 
which can be issued by the calculator are presented in 
Table III. 

TABLE III 

Command Codes for Monroe PL4 Plotter 
CODE COMMAND 



160 


X data to 


be 


passed 


161 


Y data to 


be 


passed 


162 


PEN UP 






163 


PEN DOWN 
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The plotter responds to the PEN UP and PEN DOWN commands 
only when the REMOTE switch is depressed (see Figure 10). 

In the manual mode of operation, when the REMOTE switch 
is not depressed, the pen responds to the commands issued 
by the operator. These commands are: 

WRITE (pen down) 

PEN UP 

MOMENTARY (spring loaded up, pen down while depressed) 
and are issued through switches on the plotter control 
panel. After the X data or the Y data instruction has been 
passed to the plotter, the calculator transfers the data 
from memory to the IX register for transmission to the 
plotter. See Table IV for further plotter characteristics. 



TABLE IV 

Monroe PL4 Diaital Plotter Characteristics 



Plotting Area 
Ballp oint Pen Width 
Resolution (X and Y axis) 
Repeatability 
Accuracy 
Static 
Dynamic 

Dimensions 
Width 
Length 
He igh t 
Weight 



10 inches x 10 inches 
0.009 inch 
1000 counts 

0.1 percent of full scale 

0.2 percent of full scale 
3 inch retrace will overlay 
trace within 0.035 inches 
between center lines 

17 inches 
21 inches 
7.75 inches 
48 pounds 



5 . Hewlett Packard Model 7100B Strip Chart Recorder 
The Hewlett Packard 7100B Strip Chart Recorder 
(Figure 11) has two independent pen drives and input modules 
which can accommodate analog input up to 100 V. The chart 
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Figure 10 Monroe PL4 Plotter Control Panel 
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Figure 11. Hewlett Packard Model 7100B Strip Chart Recorder. 



transport system has 12 selectable speeds and uses a 
standard chart roll with a 10-inch writing width* The 
pens move horizontally across the top of the chart each 
capable of a full 10-inch travel. Specifications are 
listed in Table V. 

TABLE V 

HP 7100B Strip Chart Recorder Specifications 

Response time 0.5 seconds full scale 

Chart Speeds 1.2 in/h; 0.1, 0.2, 0.5, 

'2 in/min 

0.1, 0.2, 0.5, 1,2 in/s 

Input Range Variable from ImV to 100 V 

Input Resistance 1 megohm 

Dimensions 

Width 16-3/4 in 

Height 8-11/16 in 

Depth 7-1/4 in 

Weight 20 pounds 

B. CALCULATOR PROGRAMMING AND DATA TRANSFER 
1 . Memory 

The calculator memory is divided into four pages. 
Each page is then further divided into columns with each 
containing 32 registers, numbered in octal from 00-37. 

The first page. Page 0, is user’s programmable memory and 
contains four columns. Each column can hold 256 instruc- 
tion codes giving a total capacity of 1024 instruction 
storage locations. The second page (Page 1) is data 
register memory also containing four columns for a total 
of 128 data registers used for main data storage. The 
third and fourth pages contain 16 columns each consisting 
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of working storage and read-only memory (ROM) . The ROM 
contains the permanently stored programs tthat perform the 
keyboard functions* The last two columns of the fourth 
page (Page 3) containing the working storage, and all of 
the first two pages are made up of random access memory 
(RAM) which can be written into and read from* Informa- 
tion stored in Pages 0, 2 and 3 is treated as a program 

instruction* A single instruction (or program step) 
occupies one byte (8 bits)* Eight single instructions can 
be stored one item per register* The fonmat is shown in 
Figure 12* Each byte contains two Binary' Coded Decimal 
(BCD) digits except the sign position* 

The calculator contains three programmable control 
registers, the Entry (E) Register, the data Address Regis- 
ter (DC) and the Index (IX) Register* The entry register 
is the primary data register* Inputs from the keyboard are 
put in the E register* Output data are formatted there, 
then transferred to a print buffer and at the end of any 
calculation the E register contains the result* The entry 
register is formatted the same as a data register with the 
addition of an overflow and an underflow digit. The address 
(DC) register contains 14 bits and contains the address of 
a specific location in memory. It is used with various 
macro instructions for defining the address to be used for 
data or code manipulation* 

The index (IX) register is a general purpose 8 bit 
register* It is divided into two 4 bit parts, IH (index 
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Figure 12, Data Register Format 
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high) and IL (index low). The IX is used as a buffer 
register for the input and output (I/O) of information to 
peripherals through the I/O bus. All information trans- 
fer between peripherals and the calculator occurs through 
the IX register. 

2 . I/O Lines 

The data output line from the calculator to the 
I/O bus is BYOT (byte, output, pin 25 on the connector). 
Data from the IX register is transmitted on this line on 
the falling edges of the clock pulses from T8-T15 time 
(see Figure 13) during an OPCB (output control byte, code 

314) command and during an OPDB (output data byte, code 

315) command. The OPCB command is used to transmit selec- 
tion and control data of peripheral devices while the 
OPDB command will send numeric data over the output line. 



The data input 


line 


into 
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pherals is BYIN (byte. 


input 


, pin 


26) . 


from the peripheral whe 
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( inpu 



305) command is executed. 
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ing. The signal OTCT (output control 
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transfer enable, pin 24) goes to a log 
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alculator from peri- 
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t data byte, code 
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OPDB (output data byte) 



(input data byte) 
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V/////////////A 
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V///////////A 



' V////////////A 



Logic 0 = +5V 
Logic 1 = OV 



Figure 13. Data Transfer Signals 
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command is executed, DTEN will remain at a logical 0 from 



T8-T15 time and OTCT will go to logical 1. 

3 • A/D Interface Data Transfer 

The operator determines which of the input signals 
to the A/D interface are converted and sent to the calcu- 
lator by entering the- proper code (shown in Table VI) into 
the computer program at decimal step numbers 130 and 135 
(see Appendix D) . During execution, the calculator loads 
the IX register with the code for the channel of the 
dependent (y) variable (decimal steps 129 and 130) then 
outputs a control byte (OPCB code 314) to set the multi- 
plexer to the proper channel then convert the signal. The 
next step (IPDB code 305) transfers the converted signal 
to the IX register. These data are then stored and the 
code for the channel of the independent (X) variable is 
loaded into the IX register and the process is repeated. 
After the two signals are in the calculator, they are con- 
verted from binary to BCD, scaled, and printed if desired, 
then passed to the X-Y plotter. 

TABLE VI 

A/D Interface Channel Codes 



CHANNEL CODE 

1 340 

2 140 

3 240 

4 040 

5 300 

6 100 

7 200 

8 000 
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4 • Digital Plotter Data Transfer 



The digital plotter has a slightly more compli- 
cated data transfer scheme because of its slower response 
time. The first 42 instructions (000-041) comprise the 
plotter subroutine and are used to control data input to 
this device. The main program branches into the sub- 
routine at the decimal steps 000, 004, 008 and 012. When 

a value is to be passed, rather than a pen command, the 
value is stored in the E register and is register 27 of 
the last column. For example, when transfer of Y data to 
the plotter is desired, the program branches to step 004. 
The program in Appendix E can be followed. 

The code for Y data (161) is loaded into the IX 
register than the subroutine jumps to octal step 16. The 
contents of the IX register is then copied into another 
register. The IX retains a copy, then the calculator out- 
puts a control byte (OPCB). The program then checks the 
lower 4-bits of the IX register (IL) to be not equal to 0. 
Since the IL register contains a 1, a jump to octal 
address 25 (decimal address 21) is performed. A RINS 
(read input status) is executed and the plotter returned 
a 1 to the H8 (left most) bit of IX if busy and the 

calculator repeats until not busy or the plotter returns a 
0. The IX register is then re-loaded with the initial Y 
code (161) destroyed by the last step then tested for a 
pen command. This test fails so the sign digit of the E 
register is passed to IL then transferred to the plotter 
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(OPDB) . The address register (DC) is loaded with the 
address of register 27 of the last column by the instruc- 
tions 273 and 277. This register contains the number in 
the E register; however, data transfer is easier from 
this location. IX is loaded with the contents of the 
address pointed to by DC which is the two exponent digits. 
These are transferred to the IX register then DC is de- 
cremented by one step (two digits) and the transfer is 
repeated. The plotter receives the exponent, the four 
most significant digits of the number and the sign of both. 
The plotter pen then moves to the point. 
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V. A/D INTERFACE PROBLEMS AND SQ3LUTIONS 



Initially the A/D interface had four features which 
inhibited it from being used with the rest of the data 
reduction system. The two major problems were (1) the 
logic for data transfer and (2) the loss of the most sig- 
nificant bit from the A/D converter. The two less serious 
and more easily corrected faults were (3) the voltage con- 
version range of the A/D converter and (4> the conversion 
of binary to decimal numbers. 

A. LOGIC FOR DATA TRANSFER 

The initial design of the A/D interfaoe assumed both 
DTEN and OCTC to remain logical 0*s throughout the entire 
normal time cycle instead of both going high from T8-T15 
time. It also assumed that DTEN would go high for data 
transfer and OCTC would go high for control output instead 
of missing the one time pulse train. Consequently the A/D 
interface could not interpret the calculator requests. In 
order to correct this design flow, the circuit labeled **A" 
in Figure 8 was incorporated. The two signals, OCTC and 
DTEN, are compared with an exclusive OR. If the two are 
different, the output of the exclusive OR is high. The 
original OCTC and DTEN are then "AND"ed to determine if 
either signal has been missed indicating the command has 
been executed. This addition changes the OCTC and DTEN 
signals to that expected by the rest of the interface. 
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B 



LOSS OF MOST SIGNIFICANT BIT 



After correcting the first problem, analog signals 
could be converted to binary and transferred to the cal- 
culator. However, only seven of the eight bits were 
being passed to the IX register. It was determined that 



the most 


significant 
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interface 


8-bit shift 


re 


gister. It was 


left 


"f loa 


ting" 


at the ou 


tput of the 


A/D 


converter then 


lost 


when 
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first tra 


nsfer pulse 


was 


being issued by 


the 


calcu 


lator . 



The circuit labeled "B** in Figure 8 was added to remove 
the first transfer pulse during an IPDB and allows the 
entire 8 bits to be transferred into the IX register. 

It uses a TTL/MST 7493 4-Bit Binary Counter and four **NOR" 



gates 


con 


figured to remove the T9 


time pu 


Ise 


from the 


clock 


to 


the 8-bit shift register 
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the 


transfer of 


data 


from 


the A/D converter. 
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First the A/D converter range was changed to OV to +4.98V 
and second the signal condition and amplifier unit was 
developed and built to enable the tape recorder signals 
to be amplified to a 5V maximum signal. By making these 
two changes the full capability of the A/D converter is 
utilized . 

D. BINARY TO DECIMAL CONVERSION 

The calculator data format is binary coded decimal 
(BCD) characters. Four bits are used to represent a deci- 
mal digit and a number is made up of by combining BCD 
characters. The A/D interface transfers a number to the 
IX register in binary representation. In order to be used 
by the calculator, this number has to be converted to BCD 
format. This is done by calculator software by a series 
of jump instruction beginning at decimal step 48 of the 
data reduction program. IX register bit LI (right most) 
is the most significant bit. If this bit is high, 128 is 
added to scratch pad register 0. If bit L2 is on, 64 is 
added, and so on such that if all bits are on, scratch pad 
register 0 totals 255. After checking all bits of the IX 
register, the scratch pad register is divided by 256. The 
result of this is to leave in the register the decimal 
equivalent of the Ov to IV analog signal converted by the 
A/C interface. 
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EXAMPLES OF CONVERSION; 



VOLTAGE 
0.25V 
0 . 50V 
0.6875V 

l.OV 



IN 



A/D OUTPUT 
0100 0000 
1000 0000 
1011 0000 

1111 1111 



CALCULATOR CONVERSION 
2^/2® = 64/256 = 0.25 
2^/2® = 128/256 = 0.50 
(2^+2®+2^)/2®= (128+32+16) /256 



= 0.6875 

(2’^ + 2® + 2®+. . .+2°)/2®= III- =.998 

Z D O 
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VI 



CONCLUSIONS 



The automatic data reduction s 
provide a quick and accurate means 
data from the XR-3 testcraft* It 
simulated and actual tape recorder 
accurate to greater than 1% of ful 
efficiently transferred into a fin 
tion and a printed record of the d 
units, can be easily obtained. Th 
its simplicity and ease of operati 
require days can now be completed 
The signal selector and interf 
capability of zeroing and calibrat 
tape recorder. It removes the sig 
vibration and can pass the data to 
or it can internally digitize the 
to the calculator for further proc 
The provision for expansion ha 
the system both through the modifi 
software. The signal selector and 
room for internal modification and 
requirements which may arise in th 
lator has over six columns of unus 
used to increase the data reductio 
future needs in this area. 



ystem has been shown to 
of evaluating performance 
has been tested with both 
signals and proven to be 
1 scale. Raw data can be 
ished graphical presenta- 



ata 


f c 


onve 


rte 


d to 


physic 


e s 


y St 


em 


i 


s s 


uper 


b for 


on . 


A 


j 


ob 


that u 


sed to 


in 


a m 


at 


t e 


r o 


f ho 


ur s . 


ace 


un 


it 


P 


rov 


ides 


the 


ing 


th 


e 


s i 


gna 


Is f 


rom the 


nal 


no 


i s 


e 


res 


ult i 


ng from 


a 


str 


ip 


c 


har 


t re 


corder , 


inf 


orm 


at 


io 


n f 


or t 


rans f er 


ess 


ing 


a 


nd/or 


plo 


t t ing . 


s b 


een 


d 


e s 


ign 


ed i 


nto 


cat 


ion 


o 


f 


hardwar 


e and 


CO 


ndi 


t i 


on 


er 


unit 


has 


addi t 


io 


ns 


to 


sat 


i s f y 


e f 


utu 


re 


f 


and 


the 


cal cu- 


ed 


memor 


y 


whi 


ch m 


ay be 


n p 


rog 


ram 


to 


fulf 


ill 



56 



APPENDIX A 



XR-3 DATA ACQUISITION SYSTEM 

The data acquisition system presently installed in the 
XR-3 testcraft is comprised of various sensors feeding 
through amplifiers and/or signal conditioners in order to 
prepare the signal for recording on a Pemco Model 120-B 
magnetic tape recorder (see Figure 14). 

A. SENSOR SYSTEM 

Sensors are presently installed to measure the follow- 
ing parameters: 



1. 


Port thru 


2. 


Starboard 


3. 


Forward s 


4. 


Aft seal 


5. 


Plenum ch 


6. 


Testcraft 


7. 


Water erne 


8. 


Pitch ang 


9. 


Pitch rat 


10. 


Roll angl 


11. 


Roll rate 


12. 


Yaw angle 


13 . 


Yaw rate 


14. 


Lateral a 


15. 


Longi tudi 



Thrust is measured by a Revere USPI-150A (30) balance 
bridge transducer. One is installed on each engine mount 
so as to measure the fore or aft thrust of each engine. 



Two pressure probes are installed in the plenum, one 
forward and one aft. Only one can be used at a time. 
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Balance Bridge 
Thrust Transducer 
(Starboard Thrust) 



Humphrey * s 
CF-18-0101-1 
Gyro 
Package 




Pitch angle 
Roll angle ^ 



— Yaw angle 



“Pitch rate 
'Roll rate - 



Yaw Rate 



Amplifier- 

Transducer 

Package 



Height 




Height 


Sensor 




Sensor 






Conditioner 



Y 1r 



±_Y-..-Y- 



Pemco 

Magnetic 

Tape 

Recorder 



Lateral acceleration Presently 

Longitudinal acceleration J not used 



Figure 14. Data Acquisition System. 
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The signal from the bridge is amplified by one of ten Grant 
Model DCAB-3 amplifiers to a range between 0,0 and 1.0 volts 
corresponding to 0.0 to 500.0 pounds thrust. Calibration 
signals of 0.0, 0.5 and 1.0 volts are also available to 

each of the amplifiers which are found in an amplifier- 
transducer package. 

Pressures are measured by pneumatic probes connected to 

Varian differential pressure transducers accurate to 1%. 

Pressure signals are amplified and calibrated similar to 

2 

thrust so as to have 60 Ibf/ft correspond to 1.0 volt. 

Velocity is measured by a vane enclosed in a water probe 
which transmits a frequency modulated signal to a Potter 
velocity conditioner. The velocity conditioner converts 
the FM signal to a dc signal from 0.0 to 5.0 volts corres- 
ponding to 0.0 to 40.0 knots. This signal is then reduced 
to 0.0 to 1.0 volt by another amplifier in the amplifier- 
transducer package. Calibration is performed by operation 
over a measured course. 

A Humphreys Model CF18-0101-1 gyro package measures 
pitch angle and rate, roll angle and rate, yaw angle and 
rate, and lateral and longitudinal acceleration. Presently 
all modes except acceleration are being used. Outputs are 
between 0.0 and 1.0 volt and range is +_ 20 degrees in pitch, 
15 degrees in roll, ^ 180 degrees in yaw, and 30 degrees 
/sec in all rates. 
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Height is measured by 

V 

Model LM4001 height sensor 
forward of the bow. Range 
^ 3 feet, both with an acc 
range is presently in use. 
is converted to 0.0 to 1.0 
conditioner . 
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B. DATA RECORDER 



The conditioned signals from the various on board 
sensors are sent to and recorded on a Pemco Model 120-B 
magnetic reel-to-reel tape recorder. The recorder has 14 
channels for electronic data plus an edge track for audio 
recording. The unit can be controlled from either the 
pilot's instrument panel or from the unit itself, located 
just aft of the pilot's cockpit. The recorder is portable 
weighing only 100 pounds and measuring 9-1/2** x 17-3/4** 

X 26-1/2**. It is easily installed and removed from the 



tester aft . 

The recorder can operate at tape speeds of 60, 30, 15, 

7-1/2, 3-3/4, and 1-7/8 inches per second and uses 

standard or precision NAB 10-1/2 inch reels and 1.0 or 
1.5 mil one-inch tape. 
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APPENDIX B 



IBM 360/67 COMPUTER DATA FITTING 
A fortran program was developed for use on 



360/67 computer to aid i 


The program 


consists of 


Part 


1: 


Ini tializat 


Part 


2; 


Data point 


Part 


3 : 


De terminati 
cients and 


Part 


4: 


Error norma 


Part 


5 : 


Evaluation 


The p 


lotting subrout 



final graphical output. It is available as a 
subroutine and utilizes a Calcomp Plotter. 

The subroutine LSF provides the fitting po 
coefficients and all the written output. It i 
cation of the subroutine LSQF2 and the procedu 
algorithm for the method can be found in Ref. 



the IBM 
design . 
Hows : 

al coeffi- 
curve 

provide the 
callable 

lynomial 
s a modifi- 
res and 

2 . 



61 



no non oooooooooooooooo 






PURPOSE 

TC FIT A CURVE TO A SET OF DATA PCINtTS USING THE 
LEAST-SCUARES METHCC kITH QRTHOGCNAL POLYNCMALS 

INPUT 

1. N=NU^'BER CF POINTS TO FIT (FCRRAT 12) 

2. THE DATA UP TO 60 POINTS (FORMAT 2F10.3) 

3. THE DEGREE OF FIT (1 THRU 20 FORMAT 12) 

4. TITLE FOR THE OUTPUT AND PLOT (FlftST 48 COLUMNS OF 
2 CAROS; FIRST CARD: TITLE, SECOND: NAME) 



OIMENSI ON FPTI200) ,XPT(200) ,X1 ( 30 , 2) , FI ( 30 , 2 ) ,X2 ( 60 ) 

X ,DEL(6C) 

REAL*:8 TITLE! 10) , LABEL/' ' / , TI TLE 2 ( 12 ) ,DMAX 

REALMS X(60),F(6G),Y(60),WI(60) ,DELY(60) , E ( 2 1 ) , SB ( 2 1 ) 
DATA WI /6C<^1. ODO/, LABEL/' E RR' / , L AEE2/ ' '/ 

DATA TITLE/10*' • / , TI TLE 2/ 12*!= • •/ 

LAST=0 
K = 5 
L=6 

READ (K,91)N 

TO READ IN THE DATA POINTS 

REAO(K, 93) (X( I) ,F(I) ,I=1,N) 

DC 7 1=1,30 
Xl( I, 1) =X( I) 

FK I, 1)=F( I) 

XI (1,2) =X ( 1+3 0) 

FK 1 ,2) =F( I+3C) 

X2(I)=X(I) 

X2( 1+30 )=X( 1+30) 

7 CONTINUE 

20 READ (K,91)KM 

IF (KH-GT..25)G0 TO 30 
KM=-KM 

8 I = l 

READ (K,200)TITLE2 
DO 12 J=l,6 
TITLE! J)=TITL£2(J) 

12 CONTINUE 

9 NPTS=30 

IF (N.LT.NPTS)NPTS=N 

IF ((I. EQ.2) .AND. ( N .LT .60 ) ) NPTS =N-30 
CALL CRAW(NPTS, Xl( 1,1 ) ,F1( 1 ,1) ,I , 1 , LABE L , T I TLE2 , 
X.2r.2,l ,0,2, 2, 5, 6, 0,LAST) 

IF (N.LT. 30)1=2 
1=1 + 1 

IF (I.LT.3)G0 TO 9 
TC PLOT THE ERROR 

CALL LSF( N,KM,X ,F ,WI ,Y,DELY, B, SB, TITLE) 

DMAX=0. C 
DC 36 1=1, N 

IF (DABS! DELYd) ) .GT.DMAX) DMAX = DAES (DELY ( I) ) 

36 CONTINUE 

DO 37 1=1, N 

DEL ( I ) =OELY( I) /OMAX/10. 

37 CONTINUE 

CALL DRAW (N,X2, DEL , 2 , 0 , LABE 1, T I TLE2, 
X.2,.2,1,0,2,2,5,6,C,LAST) 



FIGURE 15. FORTRAN CURVE FITTING ROUTINE 
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TC PLOT TH.£ FUNCTION CURVE 
KN.=-KM 

XM=(X( Ni-xa n/200.0 
XPT(l) =XUJ 

FPT(1 J=PDLY(KM,XPT,B, 1) 

DC 10 I=2t200 

XPT( I ) = XP.T( I-l) +XNT 

FPT(I ) = PCLY(KM,XPT, e, I) 

10 CONTINUE 

CALL DRAW! 200, XPT,FPT ,3 ,0 , LABE2 , TI TLE2 , EX SC , 
X.2,.2,1 ,0,2,2,5,6,0,LASTJ 
GO TO 20 

90 FORMAT (1E10.4) 

91 FORMAT (12) 

92 FORMAT (13) 

93 FORMAT ( 2F10. 3) 

98 FORMAT (• CMAX=» , FIO. 4 ) 

200 FORMAT (6A8) 

30 STOP 

ENC 
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FUNCTICN POLY 
PURPOSE 

TC EVALUATE THE POLYNOMIAL 
INPUTS 

1. N=THE DEGREE OF 

2. X= THE ARRAY OF 

POLYNOMIAL 
3- B= THE ARRAY OF 
4,. I=THE SUBSCRIPT 

FUNCTION POLY(N,X,B,I ) 
R£AL=?8 6(21), DX, C 
DIMENSION X(200) 

CX=X( I) 

C=E(N+1) 

DC 5 K=1,N 
J=N-K+1 

c=e (j)+c»Dx 

CONTINUE 

PCLY=C 

RETURN 

END 



USING NESTING 

THE POLYNOMIAL FIT 

POINTS AT WHICH TO EVAL THE 

COEFFICIENTS 
OF THE X 



FIGURE 15,(CGNTINUEC) 
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SUERGUTINE LSF 
PUPPCSE 

70 SOLVE FOR THE LEAST-SQUARES POLYNOMIAL 
INPUTS 

1. P=NUMB£fi OF DATA PCINTS(L£SS TFEN 100 J 

2. KM=TFE DEGREE OF FIT 

3. X=ARRAY OF X VALVES 
F2=ARRAY OF Y VALUES 

5. WI=ARRAY OF WEIGHS FOR THE DATA POINTS 

6. Y=GUTPUT ARRAY CF THE FIT DATA 

7. DELY=CUTPUT ARRAY OF ERROR (F2-Y) 

8. B=OUTPUT ARkAR OF THE COEFFICIENTS OF THE 
PCLYNCMIAL FIT 

9. SB=OUTPUT ARRAY OF THE ESTIMATED ERkCR IN 6 
10- T ITLE=HEACING FOR THE OUTPUT 

NOTE: X,F2tWI,Y,DELY,B»SB,AND TITLE ARE REAL=«'8 



SUBROUTINE LSF(M,KM,X,F2tWI,Y,DELY,BtSB,TITLE) 

REAL’s'S X(l),F2(l}»WI(l),Y(l)»DELY(l),b( ll,S8(Ut 
X FM tFMRtPXF tPXP,XPXP,XPXPK,PPXFP ,XP, ALPHA, BETA, PPXFt 
X F(IOO) ,P( 100),PMl ICO) tS( 100) ,ST( ICO) ,A(21,211 ,T( IOC) , 
X SJG2,SIG3,SUM£V2, FLEV , F23AR , SQ , FMF , FMKF , AM rCHI , W ( ICC) 
X , TITLE (1) ,FBAR,Xe-AR 
CC 1 1=1,21 
DC 1 J=l,21 
1 Ad ,J) = C.O 
A(l,l)=1.0 
A(2,2)=1.0 
PBAR=0. C 
SUMEV2=0. DO 
FM=0-0 
DC 5 I=1,M 

IF(WI ( I ).EC.0.0) WI(I) = 1.0DC 
5 FM=FM+WId) 

FMR=1.0D0/FM 
FBAR^O.O 
XEAR=0. 0 
DC 10 I=1,M 
K(I)=WI (I)-^'FMR 
PM(I) =OSCRT( W(I) ) 

F(I)=F2(I )>PM(I) 

SUMEV2=SUMEV2+F2( I )+WU I ) 

FBAB=F3AR + F( I )-^PM(I ) 

1C XBAR=XBAR + X( I )=*W{J ) 

T (1 ) = FBAR 
A(2,1)=-XBAR 
PXF=0.0 
PXP=0.0 

F2BAk=SUM£V2/FM 
SUMEV2=0.D0 
DC 20 1=1 

SLMEV2= SUHEV2+ W( I ) ♦ ( F2 ( I ) -F2 BA R ) 

P(I) = (X(I)-XBAR)'/<PM(I) 

PXF=PXF+P( I)vF(I) 

2C P>P = PXP+P( I )-?P( I ) 

MFM = F'' + l.C-5 
T(2)=PXF/PXP 
PMXPM=1 .GDC 
S (1 ) = PXXPM 
KMP = IA6S(KMi + i 

E(1) = T( l)-»Ad ,1)+T(2)»A(2, 1) 

E(2)=T12)>^A(2 ,2 ) 

DO 1000 K=2,KMP 
CHI=0. 

KM=K-1 

KF2=K-2 



FIGURE 15, (CONTINUED) 
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1F(KM2. EQ^O) SO TO 200 

XPXP=0. C 

XFXPM=0.0 

B(K)=0. 0 

CO 50 J = 1tM 

XP-X( J)*P( JJ 

XPXP=XPXP-»;XP*P( J) 

?0 XFXP.M=XPXPM+XP=f=PM (J) 

ALPHA=XPXP/PXP 
BETA=XPXPM/PMXPM 
PFXF=0.0 
PPXPP=0.0 
DO 100 1=1, M 
PT=P(I ) 

P(I)=( X( I J-ALPHAJ + Pdi-BETA^PMC I) 

PPXF=PPXF + P(I):«'F(I) 

PPXPP = PPXPF+P(IJ-i'P( U 
100 PP(I)=PT 

T(K)=PPXF/FPXPP 

PMXPM=PXP 

FXP=FPXPP 

A(K,1 1 =-ALPHA*A( KMl ,i )-65TA*A I KV2 , 1 ) 

A (K,KMl) = A(i<Hl,KM2 J-A( KMl.KMD^ALFHA 
A(K,KJ = 1.X» 

IF(K.LE.3) GO TO 150 
CC 12J 1=2, KM2 

12C A (K,I }=A(K/^L ,1-1 j-ALPHA=i‘A(KMl, I i - E £TA=r-A ( KM2 , I ) 

15C DO 160 1 = 1, K 

160 E (I) = B( I)+T(K}*A(K, I) 

2CC SIG3=0.0 

DC 220 I=1,M 
SC=6(K) 

D0^230^IQ=1,KKQ 
KNIC=K- IQ 

23C SC = X( n=?SQ+B(KMIG) 

YU ) = SQ 

CELYd J=Y(I)-F2(I) 

22C SIG3=SIG3-»: W( n*DELY(I )*»2 
FNKF=MFM-i<j 
SIG2=SIG3«FM/FMKF 
FLEV=( SUMEV2-SIG3J /SIG2 
EtPEV2=SIG3 
S(KJ=PXP 
CC 240 1=1, K 
24C STd) =SIG2/S( I) 

DO 300 1 = 1, K 
SE(I) = 0 .0 
DC 250 J=1,K 

250 SE( n = SB( n + ST( J }45A(J ,I)v*2 
300 SE (d = DSQRT(SE( I 1 } 

IF(KH.GT.O)GO TO 301 
CONT INUE 

IF (K. LT.KMPIGO TO 1000 
3C1 CONTINUE 

WRITE(o,95iOJ (TITLE(I),I=1, 10) 

9510 FCPiMaK IHl ,//// ,30X,10A8,// ,25X, ' COEFFICIENTS OF THE', 
X' POWER SERIES E XPA i\ S I CN ' , / / , 2 6X , ' Y (X ) =8 U ) + E (2 )^X • , 

X ' *e (3 )'->X^ f2+6{4J^X=*^3+. , . •//) 

WRiTE(6,9520) ( ( I , 8 (I ) ) , I = 1 ,K ) 

9520 FCK.XATI (15X,3( ' B I ' , 1 2 , * ) = • , 1PD12 . 5) J ) 

WFITE(6 ,9530) 

9530 FCRr^AK // ,25X,« ESTIMATES CF ERRCR FOR THE', 

X* COEFFICIENTS’,//) 

WRITEI6 ,9540) ( ( I,S8d ) ) ,1 = 1 ,K) 

954C FORMAT! ( 15X,3( ' ERRB ( ' , i 2 , ' ) =' , 1 FD9 .3 ) ) ) 

WRITE(6,9570) SIG3,FLEV 

9570 FORMAT!//, 21X, 'SUM SQ DEV = ' , IPD 1C. 3, IIX , ' F-RATIC =' 
X,1PD10.3) 



FIGURE 15„ !CONTINUED) 
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00 25 IRI=1,M 

25 CHI=CHI+0£LY( IRI i=»0ELY( IRI)/DABS( Y( IRIi) 

WRITE ,FMKF 

1111 FORMATI *C‘ ,20X, »CHISU=' f 1PD11.3 ,15Xf 'DEG CF FPEEDCP =• 
XtCPFA.O j 
WRITE (6 »9550) 

9550 FORMAT!//, 34X,' DATA FIT DIFFERENCE',/ 

2 30X,'I XII) Yin YY(I) YYI D-YI I ) ' , //) 

WRITE! 6,9560) ( ( I ,X( I) ,F2( I ) ,Y(I ) tCELYU ) ) ,1 = 1 ,M) 

9560 FCRMAT!15X,I16,3F8.3,F11.3) 
lOCC CONTINUE 

WRITE 16,9580) 

9580 FORMAT! ' 1' ) 

RETURN 

END 



F IGURE 15. (CONTINUED) 
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VELOCITY VS THRUST 5TH ORDER FIT 
COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y (X) = B( U+e(2 )*X + B(3 )*Xt*2 + B(4)=»X**3+... 



E( U= 2.57974D-03 B( 
B( 4J= 3.2S749D 02 B( 



2)= 1.10134D Cl 
5)=-5.34132D 02 



B( 3J=-9. 196960 Cl 
B( 6)= 3.26051D 02 



ESTIMATES OF ERROR FCR THE COEFFICIENTS 



ERRB( l)=8.726D-02 
ERRBi 4>=2.8160 02 



ERRB( 2)=4.646D 00 
ERRB( 51=5.5560 02 



ERRBI 31=6. 0C4D 01 
ERRBI 61=3. 913C 02 



SUM SQ DEV = 3.0510-^04 
CHISQ= 



F-RATIO = 6.944C-01 

9. 



1. 


2710-02 




DEG 


OF FREEDOM 




— DATA 




FIT 


CIFFERENCE 


I 


XIII 


Yin 


YYI I) 


YYIII-YI I 1 


1 


0.0 


0 . c 


0.003 


0. 003 


2 


0.110 


0.502 


0.467 


-0.035 


3 


0.150 


0.406 


0.452 


0.046 


4 


C.203 


0.4C2 


0.411 


C. 009 


5 


0.250 


0.411 


0.391 


-0.020 


6 


0.275 


0.403 


0.390 


-0.013 


7 


0..298 


0.398 


0.396 


— 0. 002 


8 


0-322 


0.411 


0.407 


-0.004 


9 


0-251 


0.421 


0.424 


0.003 


10 


0-372 


0.434 


0.429 


C.005 


11 


0.398 


0.443 


0.457 


0.014 


12 


0.427 


0.472 


0.477 


0.005 


13 


0.458 


0.501 


0 .499 


-0.002 


14 


0.502 


0.554 


0.539 


-0.015 


15 


C..55 0 


0.622 


0.628 


C. 006 



Figure 16. Velocity vs Thrust Output 
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VELCCITY VS THRUST 6TH ORDER FIT 
COEFFICIENTS OF THE POKER SERIES EXPANSION 



Y (XJ = 8( l)+e(2 )*X + 8(3 )=»X>»*2 + e(4 MX4 + 3+ 



B( 1)= 6,154280-04 B( 21= 1.7CC490 Cl B( 3)=-2.0364lD 02 
B( 4)= 1.CS2300 C3 B( 5-1 =-2.963280 03 B( 6)= 3.98231C 03 
B( 7)=-2.0999SD C3 8( 



ESTIMATES OF ERROR FOR THE COEFFICIENTS 



ERRB( 1)=5.957D-C2 £RRB( 2)=7.585D 00 ERRB( 3)=1.348C 02 
ERREI 4)=£.979D C2 ERRB( 5)=2.8200 03 ERR8( 6)=4.2140 03 
EPRB( 71=2. 415D 03 ERRBI 



SUM SC,C£V = 1..262D-04 F-RATIO = 7.56CC-C1 

CHISQ= 4.978D-C3 DEG CF FREEDOM = 8. 



DATA 

I X(I) Y(I) 



FIT CIFFERcNCe 
YY(I) YY(I>-Y(n 



1 O.C 0.0 0.001 

2 0.110 0.502 0.486 

3 0.150 0.406 0.434 

4 0.203 C.4C2 0.392 

5 0.250 0.411 0.392 

6 0.275 C.4C3 0.401 

7 0.298 0.398 0.4C9 

8 0.322 0.411 0.418 

9 0.351 0.421 0.427 

1C 0.272 0.434 0.425 

11 0.398 0.443 0.444 

12 0.427 0.472 0.463 

13 0.458 0.501 0.496 

14 0.502 0.554 0.564 

15 0.550 0.622 0.619 



0.001 

-0.014 

0.028 

- 0.010 

-0.019 

^ 0.002 

C.Oll 

0.007 

0.006 

-C.OOl 

0.001 

-0.009 

-0.005 

0.010 

-0.003 



Figure 16o (Continued) 
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VELOCITY VS THRUST 7TH ORDER FIT 
COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(XI = B(l).fe(2 }=J'X^e(3J*X*=^2+a(4)tX>»’»3+.. . 



E( U= 1.13850D-04 B( 2)= 2.74653D Cl BC 2)=-4.453S7D 02 
B( 4i= 3.255C0D 03 B( 5)=-l. 265470 04 B( 61= 2.7138SC 04 
B( 7)=-3. 028290 C4 8)= 1.373320 C4 B( 



ESTir'ATES OF ERROR FCR THE CCEFF ICI ENTS' 



ERRB( 1)=3. 6740-02 ERRB( 21=1.1780 01 ERRB( 31=2.6340 02 
ERRB( 4J=2.3040 C3 ERRB( 5)=1.017D 04 ERRBl 6) =2.4080 04 
ERRB( 7)=2.S18D 04 ERRB( 8i=1.420D 04 ERRBl 



SUN SQ 


DEV 


= 4-2000- 


•05 


F- 


RATIO = 9.35 


CHISQ= 


l.< 


610C-C3 




CE6 


GF FREEDOM 






DATA 


FIT 


CIFFERENCE 




I 


xm 


Y( I) 


YYUi 


YY (I J-Y( I J 




1 


0.0 


0.0 


0.000 


0.000 




2 


O.llC 


0.502 


0.49B 


-0.0J4 




3 


0-150 


0.406 


0.417 


0.011 




4 


0.203 


0.402 


0.392 


-0.010 




5 


C.250 


0.411 


0.4C4 


-0. 007 




6 


0.275 


0.403 


0.408 


0.005 




7 


0,.298 


0.393 


0.410 


0.012 




8 


0-322 


0.411 


0.412 


0. COl 




9 


0^351 


0.421 


0.417 


-0.004 




10 


0.272 


0.434 


0.425 


-0.009 




11 


0.398 


0.443 


0.444 


0.001 




12 


0.427 


0.472 


0.473 


0.001 




13 


0.458 


C. 5C1 


0.508 


0. 007 




14 


0-502 


0.554 


0.550 


-0.004 




15 


C.'550 


0.622 


0.623 


0.001 



Figure 16. (Continued) 
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VELOCITY VS THRUST 8TH ORDER FIT 
CCEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X) = B( l)+ai2)*X+B(3)*X*-<'2+B(4J*X*>i'3+.. . 



B( 1)= 8.75C39D-06 B( 2)= 4.58C880 01 B( 3)=-9.48749D 02 
B{ 4)= 8.817460 03 61 5)=-4.496C70 04 6( 6)= 1.343410 05 
B( 7)=-2. 347170 05 B( 8)= 2.222230 05 B( 9)=-8. 809140 04 



ESIIPATES OF ERROR FOR THE CCEFFICIENTS 



EPRB( 1)=1. 9730-02 ERRB( 2)=1.778D 01 ERRB( 3)=4.7730 02 
ERRBI 4)=S.1880 03 ERRBI 51=2.9760 04 ERRBI 6J=9.795D 04 
ERRfii 7)=1.£58D C5 ERRBI 8)=1.890D 05 ERRBI 9)=7.979D 04 



SUM SQ DEV = 1.0380-05 F-RATIO = 1.2190 00 

ChISC= 3.7350-04 OEG OF FREEDOM = 6. 



-. DATA 

I ( I ) Y ( I ) 



1 


0.0 


0 .0 


2 


0-1 10 


0.502 


3 


0.150 


0.406 


4 


0.203 


0.402 


5 


0.250 


0. 411 


6 


0-275 


0.403 


7 


0.298 


0.393 


8 


0.222 


0.411 


9 


0.351 


0.421 


1C 


0.372 


0.434 


11 


C-398 


0.443 


12 


0.427 


0.472 


13 


0.458 


0.501 


14 


0.502 


0.554 


15 


0.550 


0.622 



FIT 


DIFFERENCE 


YY(I) 


YYI IJ-Y( I J 


0.000 


0.000 


0.502 


-0.000 


0.407 


0. 001 


0.401 


-0.001 


0.408 


-0.003 


0.406 


0.003 


0.404 


0.006 


0.407 


-0.004 


0.418 


-0.003 


0.431 


-0.003 


0.450 


0.007 


0 .472 


0.000 


0.496 


-0.003 


0.555 


0-001 


0 .622 


-0.000 



Figure 16, (Continued) 
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VELCCITY VS THRUST 9TH CRDER FIT 
COEFFICIENTS OF THE POWER SERIES EXPANSION 
Y(X)=B( 



B( 


l) = 


9.61578D- 


•C8 


B( 


2)= 6.40666D 


01 


B( 


3)=-l .527640 


03 


B( 


4J = 


1.64A07C 


C4 


B( 


5)=-9. 964990 


04 


Bl 


6)= 3.691510 


05 


B( 


7) = 


-8.55135D 


05 


0( 


8)= 1.211070 


06 


B( 


SJ=^9.6C2140 


05 


B(10^ = 


3.2EE33D 


C5 


B( 













ESIIPATES CF ERROR FOR THE CCEFFICIENTS 



ERRB( 1)=1.807D- 


02 ERRB( 


2)=5. 


0870 01 


ERRB( 3)=1.59CC 03 


ERRB( 41=2.0680 


04 ERRB( 


5)=1. 


4690 05 


ERRE( 6)=6.263C 05 


ERRBI 7)=1.6470 


06 ERRB( 


8) = 2. 


6160 06 


ERRB( 9i=2.3C3D 06 


ERR6(1C)=8.6280 


05 ERRB( 








SUP SC DEV = 


7-2550-06 


F- 


RATIO = 1.425C-C1 


CHISQ= 2.6100-04 




DEG 


CF FREcOCM = 5. 




OATA 


FIT 


CIFFcRENCE 


I 


jcm 


Y(I) 


YY( I) 


YY(I)-Y( I ) 


1 


o.c 


O.C 


0.000 


0.000 


2 


0-110 


0.502 


0.502 


0.000 


3 


0.150 


0.406 


0.405 


-0.001 


4 


C- 203 


0.4C2 


0.404 


0.002 


5 


C.250 


0.411 


0 .406 


-0.005 


6 


0-275 


0.403 


0.404 


0.001 


7 


0.293 


0.398 


0.404 


0.006 


8 


0.322 


0.411 


0.408 


-0.003 


9 


0.351 


0.421 


0.420 


-0.001 


10 


0.372 


0.434 


0.421 


003 


11 


0.393 


0.443 


0.448 


0.005 


12 


0.427 


0.972 


0.47C 


-0.002 


13 


0.458 


0.5C1 


0.501 


-0.000 


14 


0.. 502 


0.554 


0.554 


0.000 


15 


C-550 


0.622 


0.622 


-0.000 



Figure 16., (Continued) 
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iiiifaiill III!! i|i 




VELCCITY VS THRUST lOTH CROEP FIT 
COEFFICIENTS OF THE POWER SERIES EXPANSION 
V(Xi = B( i).+ E(2]*X+8 (3)*X**2 + B(4)^X^*3 + .. . 



e( 1)=-1. 455160-07 
B( 4}=-l. 371370 C4 
B< 7)= 3.637380 C6 
B(10)=-9.S21360 06 



B( 2)= 7.792420 00 
B( 5J= 1.552C60 05 
3( 8)=-8. 564500 C6 
3(li)= 3.428700 06 



B( 3)= 4.6S988C 02 
B( 6)=-9. 642410 05 
B( 9J= 1.232150 07 
B( 



ESTIMATES OF ERROR FOR THE COEFFICIENTS 



ERRBI 1) = 1. 7770-02 
ERRB( 4)=1.CS9D C5 
ERRE( 7i=1.617D C7 
ERR6 (10=3.6710 07 



ERRB( 2)=2.0770 02 
ERRB( 5)=9.241D 05 
ERRBt 8)=3.511D 07 
ERRB( 11)=1.2280 07 



ERRB( 3J=7.2220 03 
ERRE( 6)=4.815C 06 
ERRB( 9) =4. 7620 07 
ERRE( 



SUM SC CEV = 5-6150-06 
CHISQ= i.96;0-04 



DATA 

I X(U Y(Ii 



1 


0. C 


0.0 


2 


0. lie 


0.5C2 


3 


0.150 


0.406 


4 


0.203 


0.402 


5 


0.250 


0.411 


6 


0w275 


0.403 


7 


0. 293 


0. 398 


8 


0.322 


0.411 


9 


0.351 


0.421 


10 


0-372 


C. 434 


11 


0.393 


0.443 


12 


0.427 


0.472 


13 


0.458 


C. 5C1 


14 


0.502 


0 .554 


15 


0.550 


0.622 



F-RATIO = 7.7840-02 
CEO CF FREEDOM = 4. 



FIT 


CIFFERENCE 


YY( IJ 


YY(Ii-Y( I J 


-0.000 


-0.000 


0.502 


C. CO 0 


0 .406 


-0.000 


0.4C3 


0.001 


0.408 


-0. 003 


0.404 


0.001 


0.4C2 


0.004 


0.407 


-0.004 


0.421 


—0 .Ooo 


0.433 


-0. 001 


0.448 


0.005 


0 .46tJ 


-0.004 


0.503 


0. 002 


0.554 


-0.000 


0.622 


0.000 



Figure 16, (Continued) 
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Thrust (Pounds) 



1000 




Velocity vs Thrust - 5th Order Fit 



Figure 17. Velocity vs Thrust Curves. 
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Error 



Thrust (Pounds) 



IOOOt 




Velocity vs Thrust - 6th Order Fit 



Figure 17. (Continued) 
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Error 




Velocity vs Thrust - 7th Order Fit 



Figure 17. (Continued) 
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Error 



Thrust (Pounds) 



1000 




Velocity vs Thrust - 8th Order Fit 



Figure 17. (Continued) 
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Error 



Thrust (Pounds) 



•OOOT 




Velocity vs Thrust - 9th Order Fit 



Figure 17. (Continued) 
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Error 



Thrust (Pounds) 



1000 




Velocity vs Thrust - 10th Order Fit 



Figure 17. (Continued) 
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APPENDIX C 



FLOW CHART FOR DATA REDUCTION PROGRAM 
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SYM 3 



SYM 5 
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1 

2 

•3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

:o 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 



APPENDIX D 



DATA REDUCTION PROGRAM LISTING 



nt Number 



scimal 


Code 


Symbol 


Explanation (if 


needed) 


>00 


266 




load IX register 


) 


1 


1 6 0 




with 160 


1 to transfer 


2 . 


220 




j'ump to octal | 


^ data 


3 


0 1 6 




address 16 1 


to plotter 


4 


2 6 6 




load IX with ^ 


1 


5 


1 6 1 




161 


to transfer 


6 


220 




jump to octal 


^ Y data 


7 


0 1 6 




address 16 










J 




8 


266 




load IX with 




9 


1 6 2 




162 


pen up 


1 0 


220 




jump to^ octal 


► 


1 


0 1 6 




address 16 




2 


2 6 6 




load IX with 


\ pen down 


3 


1 6 3 




163 


; 


4 


1 60 




store IX into IRQ register 


5 


3 1 4 




OPCB (output control byte) 


6 


260 




if IL of IX 7 ^ 0 


jump 


7 


02 5 


a* 


to octal address 25 


8 


3 t 4 




OPCB 




9 


220 




jump to octal 




2 0 


0 3 4 


n 


address 34 




1 


302 




RINS (read input status) 


2 


24 7 




if H8 of IX = 1, jurap to 


3 


02 5 


Ql' 


octal address 25 




4 


1 6 1 




recall IRO register to IX 


5 


3 1 4 




OPCB 




6 


24 1 




if L2 of IX = 1, jump to 


7 


05 1 


if 


octal address 51 




8 


360 




clear IX 




9 


326 




transfer digit of E to IL 


3 0 


3 1 5 




OPDB (sign transfer) 


1 


2 7 3 




load DH with 111 111 
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1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

10 

1 

2 

3 

4 

5 

6 

7 



27 7 
32 2 
3 1 5 
336 
322 
3 1 5 
336 
322 
3 1 5 

0 5 7 
066 
052 

1 1 4 
000 
06 3 
377 
24 0 
1 20 
24 1 
1 2 7 
24 2 

1 3 5 

2 4 3 

1 4 3 

2 4 4 

1 5 1 

2 4 5 

1 56 

2 4 6 
1 6 3 
24 7 
1 7 0 
1 1 1 
000 
024 
0 0 2 
00 5 
006 
0 20 
0 5 7 



er 



load DL with Oil 111 111 
load IX with (DC)* 

OPDB (transfer exponents) 
decrement DC by 1 step 
load IX with (DC) * 

OPDB (transfer first 2 digits) 
decrement DC by 1 step 
load IX with (DC)* 

OPDB (transfer next 2 digits) 



F 



RESUME (emd plot subroutine) 
syi*ibolic 52: subroutine to 






connect binary to decimal 



♦ 



0 

a 






clear dalta register 

0 

clear entry register 
no operation 
jump if LI = 1 to 
octal addres 120 
jump if L2 = 1 to 
octal address 127 
jump if L4 = 1 to 
octal address 135 



t 



0 



2 

s 

s 



jump if L8 = 1 to 
octal address 143 
jump if Hi = 1 to 
octal address 151 

jump if H2 = 1 to 
octal address 155 

jump if H4 = 1 to 
octal address 163 
jump if HS = 1 to 
octal address 170 



to reconvert 
input datta 
E register now 
^ contains input 
data valiue 

J RESUME 
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1 

2 

3 

4 

5 

6 

7 

:o 

1 

2 

■3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

lO 

1 

2 

3 

4 

5 

6 

7 



3 7 7 








3 7 7 








3 7 7 






^no operations 


3 7 7 








3 7 7 








377 








377 








3 77 






:• 


00 i 


/ 


i 


add 128 to main data register 


002 




1 


^ 0 if most significant 


0 1 0 






^ but (Ll) of IX register 


1 1 3 






is on (Ll = 1) 


000 




J 




220 




jump to octal 


062 


A 


address 62 


006 


s 




^add 64 to 14DR 0 


004 


4 




if L2 = 1 


1 1 3 








000 




«>«<' 




2 2 0 




jump to octal 


0 6 4 




address 64 


003 




1 




0 0 2. 


2 


1 


j^add 32 to IIDR 0 


1 1 3 


4 




if L4 = 1 


000 


0 






220 




jump to octal 


066 




address 66 


00 1 


/ 






006 






^add 16 to MDR 0 


1 1 3 


•f 




if L8 = 1 


000 








220 








0 70 


S 


jump to octal 


0 \ 0 




address 70 


1 1 3 




1 


add 8 to ilDR 0 


000 


0 


J 


if III = 1 


220 




jump to octal 


0 72 




address 72 


004 


4 


■) 




1 1 3 


4 


J 


add 4 to MDR 0 



(MDR) 
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1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

'0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

:0 

1 

2 

3 

4 

5 

6 

7 



000 
220 
074 
002 
1 1 3 
000 
220 
076 
00 1 
1 1 3 
000 
220 
1 00 
3 7 7 
3 7 7 
3 7 7 
3 7 7 
266 
000 
3 1 4 
305 

1 6 3 
266 
000 
3 1 4 
30 5 
05 7 
066 
07 5 
200 
0 1 4 
200 
0 1 0 
0 5 6 

0 5 7 
066 
070 

02 3 

1 1 t 

002 



lR 



if H2 = 1 
jump to octal 

address 74 

add 2 to MDR 0 
if H4 = 1 



jump to octal 
address 76 

add 1 to MDR 0 
if H8 = 1 

jump to octal 
address 100 

no operations 

next is channel select and 

data transfer code for A/D 

converter 

load IX with 

(code for Y data channel) 

0PCB)get y value in IX 
IPDB) 

s^tore. IX in IRl register 
load IX with 

(code for x data channel) 

OPCB)get x value in IX 
IPDB) 

RESUME 

} symbolic 75 point 
plot subroutine 

branch to octal 
address 14 
branch to octal 
address 10 
halt 

RESUME 

Ysymbolic sin x data scaling 
and transfer subroutine 

►EX MDR 2 for x scale 
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1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 



0 15 2 


0 2 1 




3 


1 1 1 


t 


4 


0 1 0 


• s 


5 


020 


s 


6 


220 




7 


000 


0 


8 


066 




9 


07 1 


5 


0 160 


02 3 


X 


1 


1 1 1 


t 


2 


003 


J 


3 


02 1 




4 


1 1 1 


t 


5 


0 1 1 


9 


6 


020 


~ 


7 


220 




8 


004 


4 


9 


066 




0 17 0 


.0 0 5 


S 


1 


06 5 




2 


1 66 




3 


1 1 7 




4 


003 


J 


5 


200 




6 


0 1 0 


9 


7 


056 




8 


1 1 0 




9 


002 


9 


0 180 


060 




1 


056 




2 


1 1 0 




3 


00 3 


J 


4 


0 6 0 




5 


056 




6 


1 1 0 


* 


7 


004 


4 


8 


0 6 0 




9 


05 6 




0 19 0 


1 1 0 


♦ 


1 


00 5 


S 



add X origin 



} jump to X data entry of 
plot routine 

•N symbolic arc sine Y data scaling 
and transfer subroutine 

E X KDR 3 for Y scale 



add Y oriain 



} jump to Y data entry 
plot routine 

} symbolic 5 scale factors 
storing subroutine 
advance paper 
clear flag 

set decimal point to 
3 places 

} jump to pen up entry 
of plot routine 
halt 

enter and print x scale 
for plot 

halt 

enter and print x scale 
for printout 

halt 

enter and print Y scale 
for plot 

halt 

y enter and print Y scale 



85 



1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

:0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

:0 

1 

2 

3 

4 

5 

6 

7 



0 6 0 




0 6 5 




066 




006 


s 


200 




20 1 




200 




0 5 4 


K 


1 2 7 


Br 


0 1 6 




06 7 




0 7 4 


r 


1 2 7 


Br 


06 7 




0 7 0 


s 


1 6 5 




200 




0 5 4 


'A 


1 2 7 


br 


0 1 6 




06 7 




0 5 4 


A 


1 2 7 


Br 


06 7 




0 7 1 


§ 


1 2 7 


br 


02 3 


X 


06 7 




0 7 5 


z 


1 26 


Ju 


06 7 




0 0 6 


s 


0 6 6 




6 7 4 


r 


1 1 0 


* 


00 1 


/ 


1 1 1 


t 


0 2 3 


X 


004 


4 


2 6 6 





for printout 
advance paper 



} 

} 



symbolic 6 
main routine 

branch to octal address 201 
get X and y data 



} 




branch to octal address 54 

evaluate x data 
branch if flag is set to 



symbolic r(x print 
routine 



j 

branch to symbolic sin 
> (x scale and plot routine) 



recall IRl register to IX (Y data) 




evaluate Y data 
branch of flat is set to 
symbolic 1/x (Y print 
routine) 



branch to y scale and 






plot routine 



jump to symbolic 6 again 
(start of main program) 



} 



symbolic r 
X print routine 



save X value 



scale X 
load IX with 
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1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

,0 

1 

2 

3 

4 

5 

6 

7 



3 5 4 

1 50 
1 1 1 
00 1 

0 5 7 
066 
05 4 

1 1 0 
0 0 1 
1 1 1 
023 

005 

2 6 6 

3 5 5 
1 50 
1 76 
1 1 1 
00 1 
0 5 7 
066 
0 1 2 

06 2 
060 
200 
000 
200 

004 

0 5 6 
00 1 
000 
060 
200 
000 
200 
004 
066 
000 
056 

1 1 0 
0 1 0 



K 



A 



/ 

o 



11 101 100 

print scaled x and x 

return x value 
RESUME 

} symbolic 1/x 
Y print routine 

same Y value 



scale Y 
load IX with 
11 101 101 

print scaled Y and Y 
print data 

return Y value 
RESUME 

} symbolic • plotter "zero" 
and scale check routine 



send 0 to plotter 



send 10 to plotter 



} 



syr.±)olic 0 

grid coordinate routine 



enter X grid coordinate 
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1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 



06 0 
200 
000 
056 
1 1 0 
0 1 1 
06 0 
200 

004 
066 
00 1 
056 
1 1 0 
000 
060 
056 
1 1 0 
00 1 
0 6 0 

05 6 
0 6 0 
1 1 0 
004 
056 
060 
1 1 0 
00 5 
066 
0 0 2 

0 5 6 
060 

1 1 0 
002 

0 5 6 
0 6 0 

1 1 0 

00 3 

1 1 1 
000 
1 1 1 





















t 

t 



0 



enter Y grid coordinate 



S 



4 

/ 



} 



symbolic 1 axis plot 
routine 



enter x maximum 



O 



enter x minimum 



/ 



enter y maximum 



4 



enter y minimum 



S 



} 



symbolic 2 
enter x scale 



2 



enter y scale 



J 

O 
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7 

•0 
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4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 



02 2 
00 1 
353 
1 1 1 
024 
0 0 2 
060 
1 1 1 

0 04 

1 1 1 
022 

00 5 
353 

1 1 1 
0 2 4 

00 3 
060 
056 
06 0 

1 1 0 
006 
056 
060 
1 \ 0 
007 
066 
00 3 



t 



t 

t 



t 











4 



S 



J 



e 




062 



A 



1 1 1 
00 1 
04 1 
1 2 7 
067 
070 
000 
1 2 7 
067 
07 1 
1 2 7 
067 



t 

t 

♦ 

Br 

S 

0 

Or 

0r 



compute and print length 
of X axis 



compute and print length 
of y axis 

enter interval between graduations 
of X 

enter interval between 
graduations of y 



synibolic 3 

axis plot routine 



transfer x value 



transfer y value 
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3 
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7 

.0 

1 

2 

3 

4 

5 

6 

7 

.0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

10 

1 

2 

3 

4 

5 

6 

7 



077 


a ' 




06 6 




symbolic Ig 


050 


h 


j routine to plot x axis 


04 3 


0 




1 27 






067 






0'7 0 


s 


• 


000 






1 2 7 
06 7 


ft' 




0 7 1 


s 




04 3 


0 




1 27 
06 7 


8r 




070 
0 1 2 


S 




00 1 


/ 




1 1 1 


t 




02 1 






0 1 1 
200 






004 


4 




04 3 


0 




1 2 7 
06 7 


Br 




070 
0 1 2 


S 




00 1 


/ 




1 1 1 


t 




022 


- 




0 t 1 






0 1 3 
200 


~ 




004 


4 




04 3 


0 




1 27 
067 


dr 




070 


s 




000 


0 




1 2 7 


Br 
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1 

2 

3 

4 

5 

6 

7 

:o 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 

0 

1 

2 

3 

4 

5 

6 

7 



067 

0 7 1 

1 1 1 
006 
04 1 
1 1 1 
000 
022 

0 4 3 
020 

1 26 
02 1 
06 7 

0 5 0 

1 26 
020 
06 7 

0 50 
200 
0 1 0 
000 

1 27 

06 7 
070 
062 
1 1 1 

00 5 

04 1 

1 2 7 
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TABLE VII 



DATA REDUCTION PROGRAM SCRATCH PAD REGISTER USAGE 



REGISTER 


USAGE 




AXES DRAWING SECTION 


PLOTTING SECTION 


0 


X maximum 


decimal conversion 


1 


X minimum 


temporary storage 


2 


X scale (units/0.1 in) 


X scale (V/in) 


3 


Y scale (units/0.1 in) 


Y scale (V/in) 


4 


Y maximum 


X scale (units/V) 


5 


Y minimum 


Y scale (units/V) 


6 


X interval 


not used 


7 


Y interval 


not used 


8 


grid origin X 


grid origin X 


9 


grid origin Y 


grid origin Y 
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